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ABSTRACT 



Foraging behaviors are constrained in part by the energy within resources, as 
well as the mechanical costs of procuring and consuming them. Molar enamel 
thickness is a morphological trait that differentiates African apes from hominins, 
being manifested most dramatically among hominins with large (megadont) mo- 
lars and hyperthick enamel (e.g. Paranthropus spp.). Among primates, thicker 
enamel correlates with the consumption of stiffer, deformation-resistant foods, 
probably because thicker enamel can better resist cracking under high compres- 
sive loads. Accordingly, plant underground storage organs (USOs) are thought 
to be a central food resource for early hominins, in part due to their accessibility 
and mechanical characteristics. However the potential effects of this putative 
foraging behavior on enamel wear has received little attention. Here we present a 
process-based model to investigate foraging constraints as a function of energetic 
demands and enamel wear among anthropoid primates. Our framework allows 
us to explore if and under what conditions stiff foods such as USOs are chosen as 
fallback, rather than preferred, resources. Moreover, we quantify the potential 
fitness benefits of extradentary processing or the direct employment of megadont 
molars on mechanically stiff USOs. Our model predictions bring consilience, or at 
least offer a plausible resolution, to the noted disparity between hominin isotopic 
and microwear dietary evidence. 

24 Subject headings: Hominin evolution, Enamel thickness, Paranthropus, Foraging 

models, Underground Storage Organs, Fallback foods 
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Introduction 



All animals must acquire and deliver food to their digestive systems. The mechanics 
28 of this process can result in the gradual wear, or senescence, of anatomical structures 
such as claws, beaks, and teeth. Such wear is detrimental to the foraging efficiency and 



30 reproductive success of a wide range of animals (Swennen et al. 1983 Raupp 1985; Juanes 



1992 Juanes and Smith 1995 Roitberg et al. 2005 King et al. 2005). For mammals, the 



32 oral processing, or chewing, of food is a necessary wear-inducing behavior (Stirling 1969 



McArthur and Sanson 1988; Skogland 1988 Perez-Barberia and Gordon 1998), and natural 



34 selection is expected to favor dental attributes that prolong chewing efficiency Accordingly, 
considerable attention has been focused on the microstructure, morphology, and functional 
36 ecology of mammalian molars, particularly the enamel. 

Enamel is a hard, mineralized tissue covering the elastic and vascularized dentin, and 



38 rooted by cementum to the jaws of most mammals (Lucas 2004). Oral comminution of food 



prior to digestion is, today, a uniquely mammalian behavior ( Lucas||2004 ), although gizzards 
and pharyngeal teeth can serve this function in birds and some teleosts, respectively, and 



some ornithischian dinosaurs did comminute food orally (Weishampel et al. 2004). Some 



mammals have ever-growing teeth, but primates replace their molar teeth only once, after 



which they face an adult life of wear and occasional catastrophic damage (Constantino 



et al. 2010). Thus, adult primates must find a balance between the opposing advantages of 



enamel preservation and the consumption of foods with different propensities for enamel 
wear. In this vein, the identification of contemporary form-function relationships between 
tooth enamel and diet have been instructive for inferring foraging behavior in the fossil 
record, and dental enamel has long informed debate in the discipline of paleoanthropology 



(Ungar and Sponheimer 2011). 



For example, among living mammals, relatively thick molar enamel is widely associated 
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with the consumption of stiff, deformation-resistant (hard) foods, and it follows that 
hominins such as the genus Paranthropus, which possessed large 'megadont' molars with 



hyperthick enamel, also consumed such foods (Kay 1981, 1985 Lucas et al. 2008&[a Vogel 
[et~aT1[2008| |Wood and Constantino||2007] |McGraw and Daegling||20T2| |McGraw et aLl[20T2 



Wood and Schroer 2012 Constantino 2013) Yet, the identity of these stiff food objects 



56 has long puzzled researchers and fueled hypotheses on the cost of dietary specialization 



(Baiter et al. 2012). More recently, isotopic data have emerged from a range of hominin 



taxa, including Paranthropus, that reveal a heavy dependence on C 4 -photosynthetic 
plants (which have tissues enriched in 13 C and include tropical grasses and sedges) or 



[Possibly animals that consumed these plants (Sponheimer and Lee-Thorp 1999 Ungar and 
Sponheimer||2011 ). Depending on the tissue, C4 plants can be highly resistant to fracture or 



deformation, mechanical attributes that are expected to induce wear or chip the enamel of 
molar teeth. Indeed, the molars of Paranthropus boisei are often heavily worn and deeply 



chipped (Constantino et al. 2010); and yet, paradoxically, the occlusal surfaces of seven 



well-preserved molars evince a microwear fabric that indicates a diet of soft, pliable foods 



ee (Ungar et al. 2008, 2010). These discrepant lines of evidence - indicating a diet of C4 



foods that are simultaneously hard and soft - have been challenging to reconcile, and it is 
es sometimes referred to as the C 4 conundrum. 

For Paranthropus, the consumption of 13 C-enriched tissues might have included 
70 graminivorous insects and/or the leaves, seeds, and underground storage organs (USOs) of 



grasses and sedges (Sponheimer and Lee-Thorp 2003 Sponheimer et al. 2005 Yeakel et al 



20071 |Cerling et al.|20lIHLee-Thorp|20lIt . Recently, the USOs - bulbs, corms, and rhizomes 



of grasses and sedges have attracted special attention (Dominy 2012 Lee-Thorp et al. 



2012) because they are widespread in many savanna habitats and a central food resource 



for some populations of baboons and humans (Post 1982 Barton et al. 1993 Youngblood 



76 2004 Alberts et al. 2005). Corms in particular are stiff and deformation-resistant (Dominy 
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et al. 2008 ) and a significant cause of tooth wear among savanna baboons ( Galbany et ai. 



2011). Yet the mechanical and nutritional properties of these potential foods, as well as 



the anatomical constraints of hominin dentition, are seldom factored into interpretations of 
so hominin foraging behavior, and the diet of Paranthropus remains obscured by disquieting 

discrepancies ( Grine et al.||2012 ). Here we attempt to bring consilience to these discrepancies 
82 with a model-based approach. 

The physiological and behavioral processes that yield inconsistent interpretations 
84 of diet arise can be explored with foraging models that depend on the anatomical and 
energetic states of potential foragers. Here we present a Stochastic Dynamic Programming 



(SDP) model (Mangel and Clark 1988; Mangel and Ludwig 1992 Houston and McNamara 



1999 Clark and Mangel 2000) to quantify the optimal foraging decisions for organisms that 



88 must balance energetic gain with enamel wear, while accounting for the stochastic effects of 
a variable environment. We establish our model on measurements for anthropoid primates 

90 and focus specifically on decisions affecting hominins in savanna-woodland environments. 
We determine decision matrices in which specific food resources are chosen to maximize an 

92 animal's fitness, conditional on two internal states: the organism's energetic reserves and 
molar enamel volume. 

94 This model-based approach is well-suited to test a variety of important questions about 

the effects of dental enamel on foraging, and we focus on three potentially informative 
96 lines of inquiry. First, and most essential, how does the quantity of enamel influence the 

foraging strategies of savanna-woodland anthropoid primates, and how does megadont 
98 dentition alter these strategies? Second, to what degree do these foraging decisions depend 

on resource quality and quantity, where the quality and quantity of particular food items 
ioo can vary depending on the environmental conditions? Third, can extradentary mechanical 

advantages, such as peeling, pounding, or cooking alter the influence of dental enamel, and 
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io2 to what extent do these alterations provide fitness benefits? Finally, we relate our model 

predictions to paleontological evidence of hominin diets, and conclude by showing that the 
104 model framework presented here can be used to both predict and update paleodietary data. 



Methods 

we Model Architecture 

State variables: We model the foraging decisions of an organism as a function of two 
los principle state variables: 1) the organism's net energy reserves at time t, X(t) = x, and 2) 

the organism's enamel volume at time t, V(t) = v, where time is measured in days. The 
no net energy reserves of an organism determine the amount of energy available to it for basal 

metabolic costs and daily energetic expenditures (e.g. foraging). We model a single unit of 
112 energy as 10 Megajoules [MJ], equivalent to 2388 kcal and roughly equal to the energy in 1 



kg of animal tissue (Wolfram Research 2012). Accordingly, the maximum potential energy 



lw reserves for an organism, x max , is its body size, such that a; max = 70 for a 70 kg organism, 

and a change in x by a single value is a change of 10 MJ. Because an organism cannot use 
lie all of the energy contained in its mass, we set a critical threshold for an organism's energy 

reserves (x cr it = (3/4)x max ) below which it cannot survive. By contrast, a unit of molar 
us enamel volume v corresponds to 100 mm 3 (approximately equivalent to a sphere with a 2.9 

mm radius). Among anthropoids, specific properties of molar anatomy correlate strongly 
120 with body size, and these properties can be used to approximate maximal (i.e. unworn) 

molar enamel volume, w max as a function of x max (see appendix I, fig. SI). 

122 Mechanical properties: The potential energy gained from food, as well as its impact 
on an organism's enamel, chang function of food mechanical properties. Here we 

124 consider an approximating measurement for the mechanical properties of food that takes 
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into account both the elasticity (Young's modulus, [MPa]) and the fracture toughness 
126 (Ri, [Jm -2 ]) of food i, and is generally assumed to approximate 'hardness', measured as 



y/EiRi (Lucas et al. 20086). Moreover, the toughness of a food is proportional to the 



128 amount of fiber (Lucas et al. 2000). We let r\i denote the digestibility of food % ranging 



between rji = (indigestible) and rji = 1 (completely digestible). 

Gains: For food type i, an organism's energetic gains are denoted as 7^ (in units 
of x). We identify four general food groups: 1) a nutritious, mechanically pliable, 
patchily distributed food (e.g. fruit), 2) a non- nutritious, mechanically hard, widely 
distributed food (e.g. grass leaves), 3) a moderately nutritious, mechanically hard, widely 
distributed food (e.g. USOs), and 4) a highly-nutritious, potentially hard, patchily 
distributed food (e.g. arthropods or more generally small quantities of animal tissue). 
We set the food energy density to be 717, 150, 785, and 1518 kcal/kg for fruit, grass 



leaves, USOs, and arthropods/animal tissue (Wolfram Research 2012). Because larger 



animals gain relatively more calories per foraging bout, energetic gain is calculated as 
7i = (energy density/2388) • (x max /10), where the constant (1/2388) normalizes the energy 
density of foods to units of x, and the modifier (x max /10) ensures that gain scales weakly 
with increasing body size. Although we acknowledge that macromolecules such as protein 
are vital for growth and differentiable from basic energetic requirements, here we assume 



that foraging behavior is primarily governed by caloric, or energetic, limitations (Rothman 
etalpOlU ). 



Costs: We model the daily cost of foraging for food type i, q (in units of x), as 
a function of an organism's body size, and the aggregation of food on the landscape. 
Foraging costs increase nonlinearly with body size, such that larger organisms spend a 
smaller fraction of energy per unit time than do smaller organisms. We modified the 



allometric relationship proposed by Leonard and Robertson (1997) to model daily energetic 
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150 cost (kcal/day), such that Cj = C\ ■ RMR • (l/£j)/2388, and Resting Metabolic Rate is 
RMR = 69. la^max 1 ' where C\ is the activity constant (C\ = 3.80 for moderate activity), the 

152 constant (1/2388) operates as before, and is the mean encounter rate for food i, such 
that (1/&) is proportional to foraging time. Therefore, foods that are encountered more 

154 frequently (high £j) have lower per encounter foraging costs. In addition, we assessed a 
second, costlier, version of the model, where q = {C\ ■ RMR • (1/&) + C% ■ RMR)/2388, 



156 n where Ci = 1.2, accounting for additional daily costs independent of food choice (Leonard 



and Robertson 1997). 



158 Mechanical Advantage: Extradentary mechanical advantages are introduced when an 
organism modifies the mechanical properties of a particular food resource. If an organism 

wo is incapable of altering the mechanical properties of food other than by chewing, we model 
the organism as having no mechanical advantage. Accordingly, values for toughness [Jm~ 2 ] 

162 are R = (561, 300, 265, 1345), and Young's modulus [MPa] are E = (1, 11, 5, 200) for fruits, 
grass leaves, USOs, and arthropods, respectively (Lucas 2004 Williams et al. 2005 |Dominy 



et al. 2008). We note that the high value of R for arthropods is due to their fracture 



resistant exoskeletons. 



Many animals, and primates in particular, can modify the mechanical properties of 



foods before they are consumed (Altmann 2009). We consider three potential extradentary 



168 mechanical advantages: 1) the peeling, pounding, or cooking of USOs, (i?uso = 138 and 
-^uso = 5; Dominy et al. 2008), 2) the peeling of arthropod exoskeletons (-R ar thropod = 306 



170 and -^arthropod = 0.22) (alternatively, this mechanical advantage can account for swallowing 

arthropods with minimal chewing), and 3) a combination of mechanically altering both 
172 USOs and arthopods. 

Habitat Quality: The quality of a habitat at a given time, Q(t), is variable, and this 
174 affects both the nutritional gains and foraging costs of foods. We model habitat quality 
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to be rich (Q(t) = r) or poor (Q(t) = p) at time t, such that habitat quality Q(t) changes 
176 through time in accordance to a transition probability matrix (p) where 



Prr Prp 
K Ppr Ppp / 

and, for example, p rp is the probability of transitioning from a rich quality habitat at time 
178 t to a poor quality habitat at time t + 1. The differences between rich and poor quality 

habitats are defined by differences in energetic gain, the mean encounter rate, and the 
i8o dispersion of different foods. We set energetic gain to decrease by 10% in poor quality 

habitats relative to rich quality habitats. Moreover, the mean encounter rate as well as 
182 the dispersion of food (z/j) are modified by Q(t), such that food resources are more easily 

found (higher and are less patchily distributed (higher Ui) in rich quality habitats (see 
184 appendix II for a detailed derivation of dispersion and encounter rates of foods). USOs 

are stored underground and have evolved to retain high nutrient loads during periods of 



we environmental stress (Copeland 2004). We incorporate this quality by holding the energetic 
gain, encounter rate, and dispersion of USOs constant, irrespective of habitat quality. 

188 Accordingly, habitat quality influences both food quality (in terms of its energetic gain), 
as well as the risk associated with finding food (in terms of a its mean encounter rate and 

190 dispersion). 

With this basic framework, we can asses the influence of 'wet', 'dry', and 'autocorrelated' 
192 environmental conditions on foraging decisions. Wet environments have high values of 

Prr, Ppr, and low values of p rp , p pp (such that habitat quality is generally rich). Dry 
194 environments have low values of p rr , p pr , and high values of p rp , p pp (such that habitat 

quality is generally poor). Autocorrelated environments are likely to have habitats that 
we stay rich (if they are rich) or stay poor (if they are poor), and thus have high values of 

Prr, Ppp, and low values of p rp , p pr (see table 1). We recognize that natural systems are 



198 more idiosyncratic, however this simplification allows us to assess the effects of changing 
habitat quality over time with minimal added complexity. 



200 State Dynamics 

The random variable K represents the number of food items found in a single foraging 
202 period and that with probability /W(fc)j an individual finds k items of food type i. We 

maintain this notation, upper case for random variables and lower case for specific values, 
204 for all stochastic variables. Because an organism's daily consumption is limited by gut 

volume, daily caloric gain is bounded by x s = (1/5) • x max (proportional to average 



206 anthropoid % gut volume; Milton 1989). This bound functions only to limit the energetic 
gain in the event a large amount of food (high k) is found during a foraging period. Thus, 
208 an organism's energy reserves change according to the difference equation 



X(t + 1) = X(t) + mm(kr]ai, x s ) - Q. 



(2) 



Enamel volume steadily decreases as an animal consumes resources. Although the 



underlying mechanisms of enamel loss are poorly understood (Lucas et al. 2008a), siliceous 



particulate matter is probably the most significant cause of abrasion (Lucas et al. 2012). 



212 Here we assume that hard foods (high y/E^Rl values) promote increased use of the dentition 



(cf. Organ et al. 2011), and that such use induces wear regardless of the specific cause. 



2w Accordingly, we set enamel wear, Av(Q), to be a function of: 1) the mechanical properties 
of food i and 2) a stochastic decrease in enamel volume (determined by fi). Because enamel 



216 is a nonrenewable resource, this wear cannot be undone. Teaford and Oyen (1989) showed 



that the consumption time for vervet monkeys (Chlorocebus) that ate a diet of raw Purina 
2i8 monkey chow was 8x greater than that for vervets fed on pre-mashed monkey chow. With 
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respect to enamel wear, this is equivalent to chewing 8x as much food. Teaford and Oyen 
also showed that the enamel thickness decreased by ca. 0.58 /im day -1 when vervets fed on 
the raw diet, versus ca. 0.24 /im day -1 when they fed on the pre-mashed diet. We assumed 
a linear relationship between the loss of enamel thickness (from Teaford and Oye"n| 1989), 
and consumption time, or, alternatively, the amount of food consumed, k (with a slope 
b = 0.0425). The lower-bound of this relationship (p = 0.24 /im) represents the expected 
basal enamel wear that occurs irrespective of consumption, and is used to parameterize 
the stochastic variable Q. Accordingly, given that A is the molar enamel surface area and 
Emc and Rmc are scaling constants denoting the average Young's modulus (50.44 MPa) 



and fracture toughness (1030.55 Jm 2 ) of monkey chow, respectively (Williams et al. 2005), 
enamel volume changes according to the difference equation 



250 V V EmcRmc 

v v 

Av(Q) 



(3) 



The constant 1/250 scales tooth wear to ensure the organism attains its expected longevity 



(Lindstedt and Calder III 1981), and accounts for 1) overestimation of molar enamel area 



232 (our allometric estimation includes the lateral aspects of molar surfaces), and 2) the notion 
that wear is a complex action affecting a small fraction of the occlusal surface at a given 



time (Lucas 2004). 



The basal loss of enamel thickness has an expected value E(Q) = u = 0.24 /im. As 
236 such, chewing and the daily wear unassociated with chewing ( Lucas|[2004 ) have variable 
effects on enamel wear. Specifically, enamel wear is typically small, but occasionally large, 



238 and is realized when the organism chips or looses a tooth or part of a tooth (cf. Boccara 



2004). We model the probability that Q falls within the small interval u and u + do; as a 



lognormal distribution, such that 
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fn{v) — LogNormal(c<;|u;, a) 



i exp f_o^>!, (4) 



where E(0) = ui and Var(^) = a 2 (see table 1). We assume that an individual dies when 
242 its energy reserves fall below x cr j t or its enamel volume falls below t> cr it- 



The Fitness Function 

We consider a nonbreeding interval of length T during which only foraging decisions 
influence fitness. This interval will ultimately become sufficiently large (see below) so that 



246 „we can use decisions that are independent of time (cf. Mangel and Clark 1988 Clark and 



Mangel 2000), and assume that at the end of this interval, the fitness of an individual with 



248 energy reserves x and enamel volume v is $>(x,v). For numerical computations we use 



<§>(x,v) = - ( 2 — — ) , where 

2 V x v 




v) — 0, otherwise. (5) 



The maximum fitness at t = T is realized by an organism with X(T) = x max and 
250 V(T) = fmax, and the rate of fitness decline increases as x and v approach x c „ t and v cr a. 

We explored alternatives such as <&(x,v) = (1 — ^r L )(l — ^ it ) and they had little effect on 
252 the qualitative predictions. We scaled the terminal fitness function to be 1, so it is easiest 

to consider it as survival after T for an individual whose end state is X(T) = x, V(T) = v. 

254 As described above, we also consider alternative environmental conditions (wet, dry, 
and autocorrelated), which differ in the probabilities of transitioning between rich and poor 
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256 quality conditions. We assume that natural selection has acted on behavioral decisions 
concerning diet (food choice) conditioned on physiological and environmental states. We 
258 define fitness functions 



F T (x,v,t) = max E{$(X(T),V(T))\X(t) = x,V(t) = v,Q(t) = r} , (6a) 

i 

FJx,v,t) = max E{$(X(T),V(T))\X(t) = x,V(t) = v,Q(t) = p} , (6b) 

i 

where the maximization over % chooses the food that maximizes fitness given energy 
260 reserves, enamel volume, and environment. By definition, at time T 

F t (x,v,T) = F p (x,v,T) = $(x,v). 

Maximizing Fitness by Food Choice 

262 For time periods before the terminal time t = T, an organism must survive mortality 

independent of starvation or enamel loss and choose the fitness maximizing food, given 
264 the stochasticity in food encounter. If the probability of death in a single period is set 



to m ~ e , (estimated for a subadult male chimpanzee; cf. Bronikowski et al. 2011), 



266 then F T (x,v,T) and F p (x,v,T) satisfy the equations of Stochastic Dynamic Programming 
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(SDP), such that 



F T (x,v,t) 



max (1 - m) ^ f(k T )i (prrEn^^r + minf^T^, x s ) r - (cj) r , 



Av(fi),t + 1)} 



fe=0 



+ p rp E n {F p (x T + mm(kriai,x s ) r - (ci) v ,v - Av(Cl),t + 1)} 




(7a) 



F p (x,v,t) 



max (1 - m) /(fcp)i (p P rEfi{F r (x p + min^^, x s ) p - (cj) p , 



Av(ft),t + 1)} 



fc=0 



+ PppEn{Fp(a;p + mm(kr]ai,x s ) p - (ci) p ,v - Av(Q),t + 1)} 




(7b) 



268 where the expectation is with respect to the random variable fl (eq. [3]). These equations 

identify the food i that maximizes fitness for a given energetic reserves X(t) = x, enamel 
270 volume V(t) = v, and habitat quality Q(t) at time t. Thus, optimal food selection takes 

into account the probabilistic quantities of finding K items of food i, (from the frequency 
272 distribution fx{k)i\ see appendix II) and losing a basal amount of enamel volume Q, (from 

the frequency distribution fn(oj)), as well as the probability that the environment remains 
274 in its current state or changes according to the transition probability matrix p. As eqns. 

7 are solved backward in time, in addition to obtaining the values of fitness, we create 
276 decision matrices D T (x,v,t) and D p (x,v,t) characterizing the optimal choice of food in a 

rich or poor environment given that X(t) = x and V(t) = v. Thus, two decision matrices 
278 are solved for three environmental scenarios determined by the habitat quality transition 

matrix p: 1) Wet: where habitat quality is generally rich, 2) Dry: where habitat quality is 
280 generally poor, and 3) Autocorrelated: where habitat quality shifts similarly between rich 

and poor, but tends to remain in its current state. 
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The Stationary Solutions 



As t moves backwards further and further away from T the fitness maximizing decisions 
284 become independent of time and depend only upon state, which accords with the intuition 

that far from the time at which fitness is assessed, the behavior of an organism is expected 
286 to depend on its state and on the environment, but not on the current time. We used 

these stationary decisions, which we denote by D*(x,v) and D*(x,v) for further analysis. 
288 We confirmed stationarity by calculating the summed square differences between decision 

matrix solutions from t + 1 to t, such that AD(t) = x (D{x, v, t + 1) — D(x,v,t)) 2 , for 
290 t = T — 1 to t <C T and we assumed stationary decisions had been reached when AD(t) — > 

for t <C T (for an example, see fig. S2). 



292 Forward Iteration 

The stationary solutions to the dynamic programming equations, D*(x, v) and D*(x, v), 
294 identify foods that maximize fitness for a given state (x, v) in both rich quality (Q = r) and 



poor quality (Q = p) habitats. We used a forward iteration algorithm (Mangel and Clark 



296 1988 Clark and Mangel 2000) to assess the impact that these foraging decisions have on 

the expected future fitness of individuals by iteratively solving for the state dynamics of 
298 simulated foragers over time, as specified by the state dynamics in eqns. [2] and |3j At each 

time r, the n th simulated individual, with states X n (r) and V n {r) forages for the food i 
300 determined by the decision matrices D*(X n (r),V n (r)) and D*(X n (r), V n (r)), and whether 

habitat quality is rich or poor. The time variable r is used to denote forward-iterated time 
302 units experienced by simulated individuals making foraging decisions in accordance to the 

stationary decision matrices, as opposed to the time units t used to calculate stationary 
304 decision matrices. 
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To test whether and to what extent mechanical advantages conveyed fitness 
306 benefits to hominin primates, we quantified expected future fitness, F(t\D* ,Q(t)), for 
n — 1,2, ...,N = 100, 50 kg individuals, with maximal foraging costs for days r = 1 to 

308 7~max 

= 10950 (expected lifespan of 30 years) given both the stationary decision solutions 
and habitat quality. As energy reserves and/or enamel volume decrease over the lifetime of 
310 an individual, F is expected to decrease similarly. We quantified the expected future fitness 
at time r of a population, given rich or poor habitat quality, as 

1 N 

F{t\D\Q{t)) = j^J2 F * ( X "( r )> Ur)\D*, Q(r)) . (8) 

312 Characterizing the Megadont Advantage 

We explored the potential adaptive benefits of megadont dentition using two 
314 approaches. First, we compared the proportions of foods identified to maximize fitness in 

accordance to the stationary decision solutions D*(x,v) and D*(x,v). Organisms that are 
3i6 predicted to utilize a particular resource across a greater proportion of states (x, v) may 

have fitness benefits in environments where those resources are plentiful. However, although 
3i8 the percentage of foraging choices in decision matrices is an efficient summary of potential 

dietary behavior, it should not be viewed as the proportional contribution of food to an 
320 individual's diet over time, which is calculated with the forward iteration algorithm. 

To determine whether megadont dentition provided fitness benefits over the lifetime 
322 of an individual organism, we compared expected future fitness, F, for populations of 

individuals with and without megadont dental anatomy, and during both wet environments 
324 (where rich quality habitats are more likely) and dry environments (where poor quality 

habitats are more likely). 
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326 Results 

The stationary solutions to the Stochastic Dynamic Programming equations predict 
328 that energy reserves and enamel volume have large consequences for diet choice. For 

animals without mechanical advantage in rich quality habitats, foods with the energetic and 
330 mechanical properties of fruit maximize fitness for all potential states (x,v) (fig. [l]^-). In 

poor quality habitats, fruits maximize fitness only if energy reserves were high (fig. [l]B); as 
332 reserves decline the optimal resource switches from fruit to plant USOs. Plant USOs confer 

less energetic gain than fruit, but this is mitigated by reduced foraging risk due to a more 
334 even spatial distribution (high dispersion, v). As enamel volume declined with age, the 

mechanical hardness of USOs, which produced greater enamel wear, promoted an increased 
336 reliance on riskier but mechanically pliable foods such as fruit. 

Because foraging costs scaled nonlinearly with body size, optimal foraging decisions 
338 varied accordingly. For larger animals and for each environmental scenario in our model 

(wet, dry, and autocorrelated) , a poor habitat quality was strongly associated with the 
340 consumption of riskier foods with higher nutritional yields such as fruit, whereas low 

quality, ubiquitous foods such as USOs were an important supplement (fig. [2]A) . Animals 
342 with smaller body sizes tended to rely on USOs exclusively. When habitat quality was rich, 

both smaller- and larger-bodied animals switched to a diet of energy dense foods (fruit). In 
344 the absence of extradentary mechanical advantages, extremely energy dense, but relatively 

rare foods such as arthropods were not consumed by animals of any size, regardless of 
346 habitat quality. For all mechanical advantage scenarios, as body size increased the role of 

plant USOs remained constant, however arthropods (highest nutritional gain and lowest 
348 probability of encounter) were favored over fruit (fig. [2]). Thus, in both rich and poor 

quality habitats, large-bodied anthropoid primates increased the percentage of extremely 
350 risky foods, such as arthropods, if their mechanical properties could be altered to preserve 
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enamel (fig. S3). Smaller-bodied animals lacked the energetic reserves required to forage on 
352 rare, but energy dense foods such as arthropods, regardless of their mechanical advantages. 

In poor quality habitats, stationary decision matrices reveal that megadont hominins 
354 could maximize fitness by incorporating a relatively greater proportion of fracture resistant 

foods. With no mechanical advantage, megadont decision matrices show a reduction in the 
356 percentage of fruit, and an increase in fracture-resistant USOs relative to non-megadonts 

(fig. [3]). As mechanical advantages are introduced, megadont decision matrices show 
358 similar percentages of each food item as those of non-megadonts with one important 

difference. Regardless of the mechanical advantage, megadont decision matrices have a 
360 greater percentage of USOs. 

For all populations, forward iterations revealed that expected future fitness decreased 
362 sharply early in life, but saturated as the population reached its expected lifespan of 30 

years (10950 days) (fig. [4]). Our results point to an important difference between the 
364 three mechanical advantage scenarios that are considered (none, arthropods, arthropods 

+ USOs; because the decision matrices for the USO mechanical advantage was nearly 
366 identical to the no mechanical advantage scenario, for simplicity we show only forward 

equation results for the latter) (fig. |4jA.,B; solid lines). Although the decision matrices for 
368 the arthropod and arthropod + USO scenarios are similar for non-megadont primates, these 

mechanical advantages appear to have large impacts on expected fitness. For both wet and 
370 dry environmental conditions, having either mechanical advantage provides large fitness 

benefits, but the difference in fitness between mechanical advantages is small, particularly 
372 when habitat quality is generally rich (wet conditions). 

The fitness benefits of megadont dentition are more obvious. Megadont populations 
374 have greater expected future fitness than those of non-megadonts - irrespective of mechanical 
advantage - and these differences are more exaggerated later in life (fig. |i]A,B; stippled 
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376 lines). The predicted fitness benefits generated by mechanical advantage are generally less 
for megadont populations than for non-megadont populations. 



Discussion 



Although model-based approaches have been used to investigate the foraging strategies 



of humans (Belovsky 1988), nonhuman primates (Boyer et al. 2006 Sayers et al. 2010), and 



even the interactions between the two (Levi et al. 2011), few have been applied to extinct 



primates (but see Dunbar 2005; Janssen et al. 2007 Griffith and Long 2010), and none 
have taken into account nonrenewable resources such as dental enamel. This omission is 
surprising given the functional significance prescribed to molar morphology, and its known 
impact on dietary behavior. Stochastic Dynamic Programming models demand the explicit 
expression of processes determining fitness, as well as sources of external and internal 



stochasticity (Mangel and Clark 1988 Clark and Mangel 2000). Our goal here is not to 



create an agent-based demographic model of anthropoid foraging behavior (cf. Griffith and 



Long 2010), but to assess directly the dependence of diet choice on enamel volume, and to 



quantify the extent that attributes such as extradentary mechanical advantage or megadont 
dentition alter predicted foraging behaviors. 



392 Body Size and the Role of Fallback Foods 

An important and counterintuitive result from our model concerns the effect of body 
394 size in poor quality habitats. We found that smaller body sizes can rely exclusively on 

USOs, whereas larger body sizes can subsist on a greater proportion of fruit in tandem with 
396 a smaller proportion of USOs. Such a result highlights the spatial and temporal ubiquity 

of USOs compared to fruit, and it suggests that the abundance and predictability of a 
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food resource matters more to smaller, predator-sensitive consumers with relatively higher 
metabolic demands. Perhaps compellingly, this result is compatible with the high degree 
of sexual size dimorphism observed in Paranthropus robustus - males are substantially 
larger than females - which suggests reduced food competition between females and intense 



402 male-male competition for access to females (Lockwood et al. 2007). Such social behavior 



is manifested most strongly when predation risk is high and food resources are common 
and widespread. Thus, for larger-bodied males, a diet premised on USOs is predicted to 
be suboptimal, but males are unlikely to abandon their monopolization of smaller-bodied 
females. As a result, males are expected to be more susceptible to predation due to poorer 
body condition; indeed, evidence suggests that males are disproportionately overrepresented 



in predator-accumulated fossil assemblages (Lockwood et al. 2007). 



Multiple lines of evidence suggest that plant USOs were important foods for early 
hominins. Plants with geophytic structures are both diverse and abundant in arid habitats 



(Pate and Dixon 1982 Vincent 1985 Proche§ et al. 2006), and modern hunter-gatherers 



412 utilize these resources extensively, particularly in marginal environments (Marlowe 2003 



Marlowe and Berbesque 2009). Associations between mole rats - known USO specialists 



jand hominins suggest that human ancestors lived in USO-abundant habitats (Laden and 



Wrangham 2005), while stable isotope analysis of both modern and fossil mole rats confirm 



4i6 that USO specialists have isotopic values similar to those of Australopithecus africanus 



and Paranthropus robustus (Yeakel et al. 2007). It has widely assumed that USOs served 



4i8 as fallback rather than preferred foods due to their lower nutritional content and relative 



availability (Schoeninger et al. 2001). The results of our model are in general agreement 



with this assumption, but show that the role of USOs as fallback foods varies - in part - as 
a function of an organism's energy reserves, enamel volume, and body size. 



In general, the consumption of USOs is predicted if enamel volume is relatively high 
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and energy reserves are relatively low (fig. nj). However, our model also predicts a tradeoff 
424 with respect to the role of USOs as fallback foods as body size is altered. Smaller-sized 

animals tend to preferentially use nutritious foods such as fruit in rich quality habitats and 
426 less nutritious but more ubiquitous foods such as USOs in poor quality habitats (fig. [2|\.)- 

Thus, as energetic reserves become more limiting, as they are for smaller organisms with 
428 relatively higher resting metabolic rates, fruit and USOs alternatively serve as preferred 

foods when habitat quality is rich and poor, respectively. Larger body size enables risky 
430 foraging on fruits even when habitat quality is poor. In this case, larger animals resort to 

USO consumption when energy reserves are low. Accordingly, USOs are relegated to a 
432 fallback status, and are consumed if the act of foraging for preferred foods incurs relatively 

greater fitness costs on the organism. Factors such as increased digestive efficiency are 
434 expected to counterbalance the benefits of risky foraging for very large animals, and the 

advantages these adaptations provide could be explored using a SDP approach. 

436 Although hominins are expected to have less complex hind-gut fermenting digestive 

, they were likely capable of mechanically altering their food prior 
438 to consumption, and others had megadont dentition that likely facilitated consumption of 
fracture resistant foods. We find that both mechanical advantages and megadont dentition 
440 strongly affect predicted foraging behavior. Physically altering the mechanical properties 
of resources increases potential dietary diversity by decreasing the fitness costs of fracture 
442 resistant foods. In general, the mechanical advantages (USOs, arthropods, and arthropods 
+ USOs) all increase the proportion of higher risk foods, particularly for large body sizes, 
444 in the predicted decision matrices D*(x,v) and D*(x,v) (fig. [2j3, C,D). Under the USO 
mechanical advantage, the percentage of USOs increased only marginally, while that of 
446 fruit declined. By comparison, the introduction of an arthropod or arthropod + USO 

mechanical advantage altered the decision matrices considerably. The increase in arthropod 
448 consumption was limited to organisms with larger body size, the result of larger animals 



physiologies (Milton|1989 
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having both energetic reserves and enamel volume to spare, factors that decrease the risk 
450 of obtaining rare or patchily distributed food that tend to be fracture resistant. 

It stands to reason that the evolution of megadont dentition, which increases the 
452 quantity of enamel volume for a given body size, may serve a similar function. Megadont 

dentition was a prominent morphological feature of robust hominins, though its fitness 
454 benefits are a source of contention. Megadont hominins ranged from 2.5 to 1.1 Ma, had 



smaller incisors and canines, and molars with hyperthick enamel (Shellis et al. 1998 Lucas 



et al. 2008a). Such extreme molar enamel thickness would have increased tolerance of 



fracture resistant foods (Kay 1981 Osborne 1981 Macho 1999) and it is hypothesized 



458 to have evolved in response to hard-object feeding, possibly on grass seeds (Jolly 1970) 



or plant USOs (Laden and Wrangham 2005; Sponheimer et al. 2005 Yeakel et al. 2007 



Dominy et al. 2008). Our model predicts that fracture resistant foods maximize fitness for 



a relatively greater proportion of internal states (x,v) for megadont hominins (fig. [3]). This 
morphological adaptation may have been competitively advantageous in stressed or poor 
quality habitats where fracture resistant foods are relatively more abundant. 



464 The Fitness Benefits of Mechanical Advantage and Megadont Dentition 

Both extradentary mechanical advantages and megadont dentition increase the 
466 nutritional benefits of fracture resistant foods and provide fitness advantages to hominins 

in poor habitat quality environments. Expected future fitness changes over an individual's 
468 lifetime, and is expected to diverge as enamel wears at different rates. Our process-based 

model is simplistic in that life-history stages are not included, however these simplifications 
470 enable a number of predictions regarding hominin foraging behavior: 1) behaviors that 

alter the mechanical properties of fracture resistant foods result in greater fitness; 2) 
472 these benefits are primarily realized in dry environments (where habitat quality is more 
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likely to be poor), where fracture resistant foods are relatively more abundant; 3) because 
megadont dentition results in relatively slower wear rates, megadonts have greater expected 
fitness than non-megadonts, but these benefits are primarily realized later in life. Our 
understanding of hominin diets is informed by many lines of evidence including dental 
morphology, molar microwear, and the stable isotope ratios of tooth enamel. We have 
shown that different foraging choices are predicted to maximize fitness among non-megadont 
and megadont hominins, and that these foraging strategies have different expected lifetime 
fitness values. We next show how a forward iteration approach can be used to examine 
the isotopic differences observed among hominin species, and whether the mechanical and 
physiological constraints imposed by our model are predictive of the isotopic patterns 
observed in the fossil record. 



Plants that use C3- and C4-photosynthetic pathways have different ratios of stable 
carbon isotopes ( 13 C/ 12 C). These isotopic ratios are recorded in consumer tissues and can 
be used to quantify the proportions of isotopically distinct resources in a consumer's diet. 
Animals that consume C3 plants (e.g. shrubs and trees) tend to be more depleted in 13 C 
(low S 13 C values; where R = 13 C/ 12 C and S 13 C = 1000(( J R samp i c / J R stand ard) - 1), with units of 
per-mil, %o), while those consuming C 4 plants (e.g. tropical grasses and sedges) tend to be 



less depleted in 13 C (high <5 13 C values) (Koch 2007). Here and henceforth, we refer to the 
5 13 C values of carbonate, the mineral component of bone and tooth enamel. Megadont (P. 
robustus and P. boisei), as well as some non-megadont hominins (A. africanus) , have tissues 



with relatively high 5 13 C values (Sponheimer and Lee-Thorp 1999 Ungar and Sponheimer 



2011), suggesting that they were consuming either large quantities of C4 plants or animals 



that ate C4 plants. These findings are at odds with evidence from dental microwear, which 



496 tend to support consumption of non-fracture resistant, pliable foods (Ungar et al. 2008 



2010 Ungar and Sponheimer 2011 ) 
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498 Using the forward equation framework presented above, if each food choice has 

an associated distribution of 5 13 C values, we can explore under what conditions the 
500 accumulated S 13 C value of an individual within a population changes over time as a function 

its energetic reserves, enamel volume, and the prevailing environmental conditions (see 
502 appendix III for details). Our results show that the <5 13 C values of a simulated population 

of iV = 100, 50 kg anthropoid foragers capable of mechanically altering both arthropods 
504 and USOs is influenced by both energetic reserves and enamel volume. In dry environments 

and where foraging costs are minimal, the mean 5 13 C value of simulated foragers remains 
506 relatively high (5 13 C avg ~ — 10.5%o; fig. |5]A.), due to a greater reliance on USOs (fig. S3). 

After day 3500, 5 13 C avg declines to — 11.2%o as the proportional contribution of USOs 
sob decreases and that of fruits increases (fig. (5^5). This highlights the increasing importance of 

foods that are less fracture resistant as enamel is worn - despite greater energetic costs - as 
510 well as the accompanying decrease in the mean <5 13 C value of a consumer population over 

its lifespan. 

512 If foraging costs are too great, low risk, fracture resistant foods are preferred despite 

greater enamel wear, resulting in a higher 5 13 C avg ~ — 8.8%o (fig. (5p). In this case, our 
5u model predicts S 13 C values equivalent with those observed for A. africanus and P. robustus 



(Ungar and Sponheimer 2011). In costlier environments (where energetic cost includes both 
5i6 foraging costs as well as daily costs independent of food choice), USOs tend to maximize 
fitness (fig. [5t)) until late in life, when the risks of foraging pliable but rare foods are 



5i8 superseded by those associated with declining enamel volume. 

These results help resolve the observed discrepancy between hominin isotopic ratios 
520 and molar microwear. Molar enamel is formed early in life ( Lucas|[2004 ) when food selection 



tends towards fracture resistant foods with high <5 13 C values. As enamel is worn, softer, 
522 less abundant foods with lower <5 13 C values maximize fitness (fig. |>K,B). Because fossilized 
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microwear is formed shortly before death (the 'last supper effect'), our model results suggest 
524 that softer, less fracture-resistant foods will have a disproportionately large influence on 

the microwear of teeth, particularly for older individuals. In both cases, simulated foragers 
526 incorporated foods in proportions that are not predicted by their relative abundance on 

the landscape, highlighting the importance of considering both mechanical and energetic 
528 constraints in addition to resource abundance when constructing prior distributions for 

isotopic studies. 

530 For example, our model does not show C 4 grass leaves to be solutions to the SDP for 
either non-megadont or megadont hominins (fig. [3]), despite the high relative abundance 

532 of this resource. This accords well with the aversion to C4 plants observed among 

savanna-dwelling chimpanzees ( Sponheimer et aT||2006 ) and some hominin species including 



Ardipithecus ramidus (White et al. 2009) and Australopithecus sediba (Henry et al. 2012). 



Despite their relative availability, C4 plants tend to be more fracture resistant (Boutton 



et al. 1978) and less nutritious (Barbehenn et al. 2004), possibly due to the presence of 



bundle sheath cells. Similar mechanisms have been cited to explain the avoidance of C4 



538 plants by herbaceous insects in grassland communities (Caswell et al. 1973; Boutton et al. 



1978 Pinder III and Kroh 1987). However, megadont hominins such as P. boisei have S C 



values ~ 0%o, which corresponds to a diet of 75-80% C 4 foods (Ungar and Sponheimer 



2011). Such a heavy dependence on C4 foods has led to speculation that P. boisei was 



542 possibly a grazing hominin ( Lee- Thorp ||2 01 1 Rabenold and Pearson||2011 ), an extraordinary 

idea that begs the question: under what circumstances could C 4 grass leaves maximize 
544 fitness? 

We found that when fracture resistant, low risk, and low calorie foods are hyper- 
546 abundant (such that the encounter rate of grass leaves, £ gra ss> is increased from 4 to 5; cf. 
table 1), they were represented in both non-megadont and megadont decision matrices 
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548 (fig. |6]A.,B), the consumption of C4 grasses was a fallback behavior in extremis, selected 
only when enamel volume was high and energy reserves were extremely low. Moreover, 

550 megadont dentition led to a greater percentage of states where grass leaves maximize fitness 
relative to non-megadonts (fig. [6^3) , and this is in accordance with the elevated 5 13 C values 

552 observed for megadonts such as P. boisei. 

Foraging behaviors are a consequence of both the mechanical and energetic costs of 
554 food and the constraints imposed by an organism's dentition. Dental enamel thickness 



is a highly conserved trait among individuals within modern human populations (Lucas 



556 et al. 20086), yet it varies considerably across hominin lineages in the fossil record. This 

variability is an evolutionary consequence of interactions between the dentition and food, 
558 and process-based models that integrate these ingredients can inform both the possible 

roles of certain foods as well as the potential fitness benefits of different dental morphologies 
560 or extradentary mechanical advantages. This may be particularly informative if multiple 

hominin species were conspecifics and competing for similar resources. In the future, the 
562 mechanics employed by similar SDP models could be used to explore possible outcomes 

of competitive interactions between hominin species as well as the potential advantages 
564 conferred by hunting and/or technological innovations such as cooking. 
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Fig. 1.— Stationary solutions to the fitness-maximizing equations F T (x,v) (rich quality 
habitat) and F p (x, v) (poor quality habitat) for a 50 kg anthropoid primate with no mechan- 
ical advantages. There are no differences between wet, dry, or autocorrelated conditions. 
Gray elements denote values of (x, v) resulting in mortality. 
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Fig. 2. — The proportional contribution of foods to the decision matrices of anthropoid pri- 
mates with body sizes ranging from 10 to 70 kg. Contributions of foods for the no mechanical 
advantage scenario (A.) the USO advantage scenario (B.), the arthropod advantage scenario 
(C), and the arthropod + USO advantage scenario (D.). Grass leaves are not found to be 
optimal foraging solutions in any decision matrix. Results are shown for the autocorrelated 
environment scenario; results for wet and dry scenarios were qualitatively similar. 
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Fig. 3. — Ternary diagram showing the proportional contribution of fruit, USOs, and arthro- 
pods to the decision matrices of both 50 kg non-megadont and megadont primates under 
each mechanical advantage scenario. Results are shown for the autocorrelated environment 
scenario; results for wet and dry scenarios were qualitatively similar. 
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Fig. 4.— Expected fitness trajectories for N = 100 50 kg non-megadont (solid) and 
megadont (stippled) primates over an estimated lifespan with varying mechanical advantages 
(none, arthropods, arthropods + USOs), during both wet and dry environmental conditions. 
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Fig. 5. — Forward simulation of the 5 13 C values (black line denotes running mean; gray 
band denotes maximum and minimum value), mean enamel volume, and the proportional 
contribution of food-items to the diets of iV = 100, 50 kg individuals foraging in a dry 
environment over an estimated lifespan. A. and B. When foraging costs are minimal, a 
dietary switch is observed to occur near day 3500, and labels I and II denote the pre- and 
post-diet switch. C. and D. The same simulation when foraging costs are elevated. In panels 
B. and D., the red circles and triangles denote the mean encounter rate for each food in rich 
quality and poor quality habitats, respectively. 
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Fig. 6. — Stationary solutions for the fitness-maximizing equations, F T (x, v) and F p (x, v), as 
a function of energy reserves x and enamel volume v for both non-megadont and megadont 
primates when grass leaves are hyper-abundant and for the arthropod + USO mechanical 
advantage scenario. Gray elements denote values of (x, v) resulting in mortality. 
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Parameter 


Interpretation 


Units 


Value(s): Rich quality 


Poor quality 


X(t) = x 


Energy reserves at time t 


10 [MJ] 


State Variable 




V(t) = v 


Enamel volume at time t 


100 [mm 3 ] 


State Variable 




K = k 


Number of food items found 


Count 


Stochastic Variable 




= LO 


Basal enamel wear 


[mm] 


Stochastic Variable 




7 


Energetic Gain 


10 [MJ] 


(1.5,0.3,1.6,3.2) 


(1.4,0.3,1.4,2.9) 


c 


Energetic Cost (minimal) 


10 [MJ] 


(0.7,0.5,0.7,2.2) 


(1.1,0.5,0.7,2.2) 




Energetic Cost (maximal) 


10 [MJ] 


(1.4,1.2,1.4,2.8) 


(1.8,1.2,1.4,2.8) 




Mean encounter rate 


time -1 


(3,4,3,1) 


(2,4,3,1) 




Dispersion 


NA 


(3,5,3,2) 


(2,4,3,1) 




Digestibility- 


NA 


(0.9,0.7,0.8,0.9) 


Same 


A 


Molar surface area 


[mm 2 ] 

L J 


^ -1 7T L r> . 


Same 


b 


Slope of enamel wear 


[mm/k] 


0.0425 


Same 


E 


Young's modulus 


[MPa] 


(1,11,5,200) 


Same 


R 


Fracture toughness 


[Jm- 2 ] 


(565,300,265,1345) 


Same 


U) 


Expected basal enamel wear 


fim 


0.24 


Same 


a 


Basal enamel wear SD 


fim 


1.6 


Same 


d 


Prob. of death at time t 


NA 


4.54 x 10~ 5 


Same 


Q(t) 


Habitat quality at time t 


NA 


r 


P 


n 


Quality transition probability 
matrix: (p rr , p rp ; p pr , p pp ) 




Wet (0.8, 0.2; 0.2, 0.8) 
Dry (0.2, 0.8; 0.8, 0.2) 
Auto. (0.8, 0.2; 0.8, 0.2) 




$ 


Terminal fitness function (f = T) 








F 


Fitness function (t < T) 








D*(x,v) 


Stationary decision matrix 








F 


Expected future fitness 









Table 1: Parameters and variables in the dynamic state variable model. Parenthetical values 
(except for p) refer to those for foods: (fruit, grass leaves, USOs, arthropods). Values for 
E and R are those when no mechanical advantage is included. See methods for relevant 
references. Auto. = Autocorrelated. 



